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Abstract. For A</< 4 models, the introduction of a large field cutoff improves 
significantly the accuracy that can be reached with perturbative series but the 
calculation of the modified coefficients remains a challenging problem. We show 
that this problem can be solved numerically, and analytically in the limits of large 
and small field cutoffs, for the ground state energy of the anharmonic oscillator. 
For the two lowest orders in A, the approximate formulas obtained in the large 
field cutoff limit extend unexpectedly far in the low field cutoff region and there 
is a significant overlap with the small field cutoff approximation. For the higher 
orders, this is not the case, however the shape of the transition between the 
small field cutoff regime and the large field cutoff regime is approximately order 
independent. 



PACS numbers: 11.15.Bt, 12.38.Cy, 31.15.Md 

1. Introduction 

Perturbative methods and Feynman diagrams have played an important role in the 
development of quantum field theory and its applications. However, perturbative series 
usually have a zero radius of convergence 0. For scalar models with A0 4 interactions, 
the coefficients of perturbative series grow factorially. For any fixed, strictly positive, 
value of A, there exists an order beyond which adding higher order terms diminishes 
the accuracy. This feature will restrict our ability to perform high precision tests of 
the standard model (for instance, g — 2 of leptons and the hadronic width of the Z°) 
during the next decades. 

The large order behavior of the series is dominated by large field configurations 
which have little effect on low energy observables. Introducing a large field cutoff 

E] m the path integral formulation of scalar field theory, dramatically improves the 
large order behavior of the perturbative series. In two non-trivial examples 3 , this 
procedure yields series that have finite radii of convergence and tend to values that 
are exponentially close to the exact ones. This also allows us to define the theory for 
negative or complex values of A, a subject that has raised a lot of interest recently 

01 Ej ■ An important feature of this approach is that for a perturbative expansion at a 
given order in A, it is possible in some cases to determine an optimal field cutoff using 
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the strong coupling expansion [H1E|, bridging the gap between the two expansions. In 
other words, the modified perturbative methods allow us to take into account non- 
perturbative effects. 

Despite these promising features, calculating the modified coefficients remains a 
challenging technical problem. While developing a new perturbative method (see e.g. 
Ref. [HI E] for the 8 expansion) , it is customary to demonstrate the advantages of a 
method with simple integrals and the non-trivial, but well-studied|l()l IT], case of the 
anharmonic oscillator. A simple integral has been discussed in Ref. (i . In this article, 
we show not only that this program can be completed in the case of the anharmonic 
oscillator, but also that the results show remarkable properties: 

• For the two lowest orders, the approximate formulas obtained in the large field 
cutoff limit extend unexpectedly far in the low field cutoff region. 

• For the higher orders, the transition between the small field cutoff regime and 
the large field cutoff regime can be approximately described in terms of a single 
function. 

The results are presented in the following way. In Sec. |3 we define the model 
considered and the numerical calculation of the modified coefficients. In Sec. we 
discuss the radius of convergence of the modified series from a numerical point of view. 
In Sees. 0Jand|31 we present the approximate methods at small and large field cutoff. 
Rigorous bounds on the radius of convergence of series encountered in the previous 
sections are given in Sec. The question of the interpolation between the small 
and large field cutoff regimes is discussed in Sec. [7| The importance of a complete 
understanding of this question, as well as the extension to field theory are discussed 
in the conclusions. 



2. The problem and its numerical solution 

In this section, we introduce the anharmonic oscillator with a "field cutoff" and 
we explain how to calculate numerically the perturbative series in the anharmonic 
coupling A. This method was first used in |3j and tested with the known results |TU) 
up to order 20. It is a perturbative version of the numerical method proposed in 
Ref. |1 1 j . For convenience, we use quantum mechanical notations instead of field 
theoretical ones <j> — ► x, m — > ui and the field cutoff will be denoted x na x- We also 
use units such that % = 1 and the "mechanical mass" to is 1. However, the harmonic 
angular frequency u), will sometimes be used as an expansion parameter and will be 
kept arbitrary in the equations. In these units, x max ^JTo and X/tu 3 are dimensionless. 
The hamiltonian reads 

H=£ + V{x) (1) 

with 

yr x ) = [ Wx 2 + \x 4 if \x\ < X max , 2 « 

\ 00 i f \ x \ > W 

Our main goal is the calculation of the modified coefficients Eq (y/ojx max ) of the 
perturbative series for the ground state energy: 

oo 

E (x max ,uj,X) =uj^2E^ k) (y / ujx max ) x (\/uJ 3 ) k , (3) 

k=0 
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For this purpose, we will solve perturbatively the time independent Schrodinger 
equation with the boundary condition &(x max ) = 0. We proceed as in Ref. 
Setting 

oc e Jj: o , (4) 

the Schrodinger equation reads 

L' + 2(V - E)K + GL = (5) 
K' + L + GK = (6) 

where G(x) is an arbitrary function. The second of these equations implies that 

y(x)ocK(x)er d y G ^ (7) 
It possible to show that if G and V are polynomials: 

• Eqs. lloKil) can be solved by power series which define entire functions; recursion 
formulas and initial conditions (different for even and odd solutions) are given in 
Ref. HU- 

• the zeroes of $ are the sames as the zeroes of K. 

This implies that the energy levels can be obtained by solving 

K(x max ,u,X,E) = (8) 

for the variable E, and that polynomial approximations can be used for this purpose. 
We now use the perturbative expansion 

oo 

E^u J Y / E^(X/u J 3 ) k . (9) 

fe=0 

We assume that u > 0. The case u> — is simpler and will be discussed in Sec. 
El By construction, the coefficients E^ are dimensionless and can only depend on 
\fux max when their value is fixed using Eq. (JSJ). In the following, this dependence 
will sometimes be kept implicit in order to reduce the size of some equations. We also 
recall that if we had not chosen the units H = m = 1, the dimensionless quantities 
used above would read (^u>m/hx max ) and hX/m 2 u 3 . 
We need to set the expansion 

K{x max ,u,X,E) = K^\^x max ,E^)+ (10) 

(X/u 3 )K^(V^x max , EM, EV>) + ... (11) 

equal to zero order by order in A. We can approximate the by polynomials. 

We then start by solving /^(v^x™^, £(»)) = for E^ using Newton's method. 
This correponds to the problem of an harmonic oscillator potential that becomes 
infinite at x = ii mnl . The various zeroes of the polynomial correspond to the 
(even or odd depending on how we construct K |11| ) spectrum of this model. The 
ground state energy is obtained by taking the lowest even solution. By increasing 
the degree of the polynomial approximation, we can stabilize the numerical value 
with great accuracy. The independence on G of the exact equations, before we 
use polynomial approximations, can be used to test the numerical accuracy of the 
polynomial approximations used for K. In practice a good choice is made by having 
K of order 1 for E close, but not fine tuned, to its actual value. Since the potential 
is even, it is natural to have K even and G and L odd, however we noticed that 
introducing a small parity breaking in G usually improves the numerical stability. We 
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Figure 1. The numerical values of i?fc(a; max ) for k = 1, . . . 10 and lu = 1. 

then solve {x max , E^ , E^ ) = for . Since E^ appears only linearly in the 
order A expansion of K, it is a linear equation for this quantity and it can be solved 
trivially. The same reasoning shows that the higher order equations are linear in E^ . 

By using this method, we have calculated the first ten coefficients for u> = 1 and 
values of x max between 0.5 and 7. In order to allow a comparison among the different 
orders, we define the ratios 

RkiV^Xmax) = E^\y/uJX max ) / E ( Q k) (oo) , (12) 

which all tend to 1 in the x max — > oo limit and to in the x max — > limit. The 
numerical values of these ratios are shown in Fig. ^ A striking feature is that the 
curves for the various orders have approximately the same shape and that we could 
approximately superpose them by appropriate horizontal translations. Before studying 
these curves in the large and small x max limits, we will discuss numerical estimates 
for the radius of convergence of the modified series, that can be extracted from this 
data. 

3. Numerical evidence for a finite radius of convergence 

If we calculate the integral of a function defined by a series which is uniformly 
convergent over the range of integration, it is legitimate to interchange the sum and 
the integral. On general grounds, one would thus expect that if we calculate the 
partition function of a properly regularized theory (say on a finite lattice) with a 
field cut, the perturbative expansion becomes convergent for arbitrary complex values 
of the coupling. The ground state energy of the anharmonic can be obtained from 
the logarithm of the partition function and consequently, we would expect that its 
perturbative expansion will have a finite radius of convergence. This radius will depend 
on the position of the complex zeroes of the partition function. 

It is sometimes possible to estimate the radius of convergence by considering the 
empirical asymptotic behavior of the perturbative series. However, with a limited 
series, transient behavior is often encountered. It is clear from Fig. ^ that if x max is 
large enough, the beginning of the series looks like its x max — > oo limit which grows at 
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Figure 2. \n\E^ (x max )\ versus k for x ma x=0.5, 0.7, ...1.3 (from bottom to 
top); u) was set to 1. 
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Figure 3. lnA c (x max ) versus \nx max for u> = 1. The line is 4.22 — 5.89x. 



a factorial rate. Consequently, in order to attempt to probe the asymptotic behavior, 
we will only consider x max such that the available coefficients are significantly smaller 
than their value when x max is infinite. For this reason, we have limited the range of 
investigation to x max < 1.3. If we have a finite radius of convergence A c , we expect 
that for k sufficiently large, 



\Eq '(\/uJX m ax)\ OC (P(\/ujx m ax)) 

When this is the case, X c (x max ,ui) — u> 3 / P(-\/uJx max ). 



(13) 

In Fig. we have plotted 

Ih\Eq K ' (y/ujx max )\ versus k. We see a clear linear behavior and the slope can be 
interpretated as lnP(^/ujx rnax ) . We can then study A c (x max ,u) as a function of x max . 
This is done in Fig. [3] where a log-log plot shows a clear linear behavior. The slope 
is near 5.9 which is close to the value 6 that, we will argue in Sec. El should hold in 
the limit of small x max (see Eq. I|4()|l). We conclude that in this limit, our numerical 
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data suggests 

A c (a; max ,cj = 1) ~ 65 X x m \ x . (14) 
4. The large ir max limit 

In this section, we first discuss the harmonic energy spectrum with vanishing boundary 
conditions at ±x max and then treat the anharmonic interactions with the usual 
perturbative methods. We work in the limit where x max is large. This means that 
if we consider the n-th energy level, x max should be much larger than the largest 
zero of H n (xy/Iu). According to Eq. (6.32.5) of Ref. J2jj this imposes the restriction 
Xmax >> y/2n/u! for large n. When all the zeroes oiH n are well within [—x max , x max ], 
the new boundary conditions require changes that are exponentially small. It is clear 
that for any fixed x max , this condition will be violated for n sufficiently large. 

4-1- The harmonic case 

We first need to calculate the energy eigenvalues at A = 0, En\y/^x max ), and their 
corresponding wave functions ^>n (x). The wave function always depends on x max , 
but this will be kept implicit. We impose the boundary conditions (0) = for 
n even, ^n\o) = for n odd, and ^^ \x max ) = in both cases. We will not 
pay attention to the overall normalization until the end of the calculation. When 
i max — > oo, we have En\\foJx max ) — (n + 1/2) and we will use 

e n (-y/ujx max ) = E^\sjujx max ) - n - 1/2 (15) 

as our expansion parameter. Again, e„ is a dimensionless quantity that can only 
depend on -Jujx m ax and we will often keep this dependence implicit. We write 

^\x)^K^\x)e— 2 / 2 . (16) 
This corresponds to a choice G = —lux in Eqs. (151611 and consequently 

- (l/2)tf(°>" + loxK^' - nuKW = e , nLjK W . (17) 

We then expand 

K^=K^+e n K^ + ... (18) 

When x max — > oo, e„ — > 0, and we want to recover the usual harmonic oscillator 
solution. Consequently, we set Kn^°\x) = H n (y/ujx), the n-th Hermite polynomial 
which is obviously a solution of Eq. (|17fl at order in e n . At order e„, we have 

- (l/2)if< >W" + uxKWW _ n^Wti) = uHn (rf5x), (19) 
with the boundary conditions 

#(°)«(0)=i4°)«'(0)=0. (20) 

Remarkably, this inhomogeneous equation can be integrated exactly in two steps. 
First, we set 

K^(x)=H n (V^x)G n (x) . (21) 
This removes the explicitly n-dependent term and the equation depend only on G'\ 
- (l/2)H n (Vux)d„+[uxH n (Vwx) - H^{V^x)Vlo] G n = cvH n (^x), (22) 
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We then write G'(x) — N(x)e LUX which removes the term linear in x and allows us to 
write the 1. h. s. as a total derivative. The solution is then 

G n {x) = -2oj [ dy(H n (V^y)y 2 e^ 2 C dze^ (H n (^z)) 2 , (23) 
Jo Jo 
One can check that the lower bounds of integration imply the boundary conditions 
of Eq. I|2U|) . This is obvious when n is even. When n is odd, the z integral is of 
order y 3 when y — > which overcompensates the y~ 2 of the other factors. Note that 
G n (x) is independent of x max - The integrand of the y-integral has a double pole 
at the zeroes of H n , however it has no single pole. This can be seen by plugging a 
Laurent expansion of G' with poles of order 1 and 2 about a zero of H n , and using 
the relation between the first and the second derivative of H n at this zero provided 
by the defining equation for Hcrmitc polynomials. This forces the coefficient of the 
simple pole to be zero. Consequently when we express G as the integral of G', we 
can go around the double pole cither above or below the real line and obtain the same 
result. In other words, we can regularize the y-integral by replacing (H n (y / ujy))~ 2 by 
(H n (y/ujy) ± ie)~ 2 . From these considerations, we see that G(x) develops a simple 
pole when x approaches a zero of H n , which compensate the simple zero of H n . Due 
to the absence of single pole in G' , there is no logarithmic singularity in G. 

e n can be calculated by imposing the condition that the wave function vanishes 
at Xmax- This translates into the simple equation: 

^-n^/^^max) — 1 G n [x max) • (^4) 

As e ra increases, the non-trivial zeroes of the wavefunction move toward the origin. 
The shift of these zeroes is approximately linear in e n for large x max - For instance, 
for n =2, the shift of the non-trivial zero X($\ obeys the approximate linear relation 
Ax (0) ~ -0.145e 2 . 

We now discuss in more detail the case of the ground state (n = 0). It should be 
noted that in this case the double integral can be expressed ^Sj m terms of generalized 
hypergeometric series 

G (x) = -lux 2 2*2(1,1; 2,3/2; lux 2 ) (25) 

2 ^ kr(k + 1/2) ' 

k— 1 

This series defines a function that converges over the entire complex plane. The 
asymptotic form of Go(x max ), for x ma x large, can be worked out by noticing that the 
integral over y in Eq. l|23|l comes mostly from the region y ~ x max and consequently, 
we can extend the integral over z to infinity with exponentially small errors. The 
integral over y can then be approximately performed by expanding the argument of 
the exponential about x max - Using Eq. I|24l) . we conclude that asymptotically 

eo{\/ujx max ) ^ 2^-x ma xe~ UJX ™^ (26) 

This result coincides exactly with the semi-classical estimate given in Ref. |TI]. The 
validity of Eq. I|2t)|l can be checked numerically. For instance, for x rnax — 7 and u> = 1, 
we get e = 4.14 x 10~ 21 from Eq. (|26*|) which agrees well with 4.10 x 10~ 21 obtained 
numerically. 

In Fig. 0J we have plotted the asymptotic form Eq. but also the full integral 
form 

eo(V^w) = V(2w / dye^ 2 f ' dz^) , (27) 
Jo Jo 
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Figure 4. Numerical values of Ro(x ma x) (crossed circles) for cj = 1, the integral 
formula of Eq. 1271 (solid line), the asymptotic Eq. 1261 (dash-dot line), the 
leading term (<x x^ ax , dash line) of the small x max expansion studied in Sec. 
151 and together with the correction (dots). For x max > 1.5, the integral 

formula, its asymptotic form and the data are undistinguishable with the eye. 



One notices that the integral form stays accurate to much lower values of x m ax, even 
at very low values of x ma x , where it has no reason to be accurate, it gives a reasonable 
answer. This observation will be discussed further in Sec. 



4-2. Anharmonic corrections 

Having solved the problem for A — at first order in e„, we can use the usual methods 
of perturbation theory. Eq can be calculated by taking the average of x* with the 
corrected ground state wave function constructed above: 

E { 1] = ^N- 1 I*™* dx\¥ °\x)\ 2 x 4 . (28) 
Jo 

with the normalization constant 

N= [ dx\y ( 0) (x)\ 2 . (29) 
Jo 

Proceeding as for Eq. I|26() . we obtain in leading order 

~ (3/4) - Tr^w^^e-^ (30) 

Note that in this calculation, we have replaced Go (x) in the integral by its asymptotic 
form (oc x~ 1 e ulx see Eqs. (|24l ^6|0. The exponentials exactly cancel inside the 
integral and we are left with the integration of x 4 /x. This is justified by the fact that 
most of the contributions come from the large x region. The validity of the Eq. (|30|l 
can be checked numerically. For instance, for x max = 7 and u> = 1, Eq. I|30|) gives a 
correction —4.97 x 10 -18 while numerically we obtain —5.01 x 10 -18 . 
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Figure 5. Numerical values of Ri(x ma x) (crossed circles) for lj = 1, the integral 
formula of Eq. 1281 (solid line), the asymptotic Eq. 1301 (dash-dot line), the small 
%max expansion studied in Sec. 151 truncated at order 4, 8, ... 24 (dots). Series 
analysis suggests that the small x max expansion converges for x max < 2.9. 



The good accuracy of Eq. (|3U|l is shown in Fig. |S] In this figure, we also show 
the full integral formula of Eq. i|28|) . where the integral was done numerically using 
the square of the order e approximation for the ground state wave function, without 
discarding order e 2 terms in the square. Again the integral formula stays valid at 
unexpectedly low values of x ma x- This will be explained in Sec. 

The derivation of Eq. (|30|) is quite simple: eo give a contribution proportional to 
imaie~ WIm " and the ^-integral gives a factor x1^ ax . We conjecture that at leading 
order, a similar situation is encountered for higher order terms, and that 

1 - i? fe (V^max) « a:£+ 1 e- W9 £- . (31) 
For this purpose we have studied numerically the quantity 

Qk{V^x max ) = e +ux ™*(l - i4(V^ma X )) , (32) 

which according to the conjecture should scale like x^^x ■ I n Fig. El we have set 
uj = 1 and displayed ln(Qk(x max ) versus hx(x max ). The solid lines represent the linear 
functions Ak + (4fc + l)ln(x max ) . The constants Ak have been fitted using the last 
ten data points. Their numerical values are Aq — 0.802 which compares well with 
the prediction of Eq. H2tjfl ln(47r -1 / 2 ) ~ 0.814, and A\ = —0.276 which compares well 
with the prediction of Eq. ln(47r" 1 / 2 /3) ~ -0.285. The linear behavior of the 

higher orders supports the conjecture. 

5. The small x max limit 

When the field cutoff x max — > 0, the potential term in the harmonic oscillator is much 
smaller then the kinetic term, and in first approximation we can use the free particle 
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in a box of length 2a; max . The energy levels diverge as x„? ax and it is convenient to 
introduce the following rescalings 



X I Xmax 



X = \r° 
" /x,Xj max 



We then obtain a new hamiltonian 
-(- 



1 , d , , 
H = —^? + Vi (33) 



with 



With these rescalings, x max has disappeared from the problem and we can now expand 
E as a double series in A and w. This can be done using the usual methods of 
perturbation theory since we can solve the problem exactly when A = a) = 0. The 
perturbative series becomes 

oo 

E n = E^X k Cj 21 , (35) 

with dimensionless coefficients E^' 1 ' . Scaling back to the original problem, we obtain 

oo 

E n = £ E^XWx*!**-* (36) 

k=0,l=0 

Comparing with Eq. 0, we conclude that with this expansion 

oo 

4 fc) (^Lj = 5><*'° x MLJ 2 '^- 1 ■ (37) 

Z=0 
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Table 1. The coefficients E\ 
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Figure 7. ln{\E (k ^\) versus i, the order of expansion in ui , for k = 0, 1, . . . 5. 
The lines represent linear fits for a given k. 



Obviously, this implies that the asymptotic behavior when x max — > 0, for the ratios 
displayed in Fig. ^ is 

Rk(V^ x max) oc x^ax ( 38 ) 
We checked numerically that this asymptotic behavior is correct. Using the data of 
Fig. one can indeed determine the constant of proportionality Eq 1 ' ' 1 calculated 
independently below with more than two significant digits. 
For the ground state, 



,(0,0) 


n 2 




o — 


8 




,(1.0) 


1 


4 


o — 


5 ~ 




,(0,1) _ 


1 


1 


o - 


6 ~ 


—> 



24 

^4 



E r' = - 6 -- 2 (39) 

Higher orders require infinite sums. In practice, since we see from Eq. 134|l that the 
calculations can be performed with x max — 1, it is easier to proceed numerically as 
in Sec. [5] but with a double series. The numerical results are given in Table ^and 
displayed in Fig. [7| 
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Figure 8. Estimated values of x ma x c{k) for which we reach the radius of 
convergence for u = 1 (empty boxes) and values of xo(k) (black circles) where the 
crossover occurs as defined in Sec. |7| for k = 1, ... 5. 



The approximately linear behavior of Iu(\Eq |) in I at fixed k shown in Fig. \7\ 
suggests that the power series in w 2 xf nax , for the coefficients of (A/w 3 ) fc have a finite 
radius of convergence (see Eq. As the slope vary from -4.7 for k = to -3.3 for 

k — 5, we expect that the series converge for \x max ^/uj\ < e 4 7 / 4 ~ 3.3 for k = to 
e 3.3/4 ^ 2.3 f or fc — 5. Fig. [S] shows a steady decrease of the radius of convergence 
when k increases. At the same time, the value of x ma x where the transition between 
the small and large cutoff limit occurs, denoted x${k) increases with k. This quantity 
is defined more precisely in Sec. [7| Consequently, it seems clear that as the order 
increases, a gap opens between the range of validity of the small x ma x expansion and 
the crossover region. 



6. Rigorous lower bounds ont radius of convergence 

In Sees. and 03 we have discussed the radius of convergence of series from a 
numerical point of view. We now return to this question with more analytical 
considerations. First of all, it is possible to get an order magnitude estimate for 
the radius of convergence by using a simple argument. If we look at a typical term 
of the perturbative series at order A fc , it is the product of k matrix elements for the 
perturbation potential x 4 divided by k — 1 differences of energy level. The matrix 
elements of x 4 are bounded in magnitude by x max . The differences of energy depend 
on x max - If x max is very small, in other words if x^^w 2 << w in units where the 
mechanical mass and Ti are set to 1, the harmonic term is a perturbation and the 
energy differences are proportional to x m 2 ax . For instance, the difference between the 
ground state and the first excited state will be approximately 3ir 2 / (8x max ) . In the 
opposite limit of large x max , the level are approximately equidistant with a separation 
u>. From this we expect that the perturbative series converges for |A| < \ c {x m ax,^>) 
with 



x 



" 6 if x << lo^ 1 / 2 



1 U)X m * x if Xmax » OJ L ' Z 
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This argument does not take into account aspects such as the number of terms 
present or their signs, however, a rigorous lower bound on the radius of convergence 
can be obtained by simply using one half of the difference between the ground state 
and the first excited state of the unperturbed system in the above argument. This is 
formulated precisely in Theorem XII. 11 in Ref. ^3] and proved in Ref. ^S]- The final 
result in our case is 

\ c (x max ,u) > (l/2)(^ 0) - E^)/x max (41) 

Note that Theorem XII. 11 is more than what we need here because the perturbation 
is a bounded operator and we do not need to compare its norm with the norm of H . 

In the case to = 0, dimensional analysis dictates \ c {x max ,uj = 0) = Cx max for 
some positive C, and Eq. (|41|l implies 

\ c (x max ,u; = 0)> (3Tr 2 /m)xx~ 6 ax (42) 

This bound can be compared with numerical estimates obtained from the series with 
coefficients E^' ^ defined in Sec. |S] and which only requires one calculation with 
Xmax = 1 instead of the two step procedure of Sec. Using the coefficients for 
k = 5, . . . 10, we obtain 

K(x max ,uJ = 0) ~ 53 X x m 6 ax , (43) 

which is quite close to result obtained with u> = 1 at low x max in Sec. The rigorous 
bound on C, 37r 2 /16 ~ 1.85 is about 30 times smaller than the empirical value 53. 
In summary, the lower bound is compatible with our numerical estimate, but is not 
sharp. 

The same argument can be applied for the calculation of £q ^ as a power series in 
(l/2)u! 2 for the perturbation (1/2)lo 2 x 2 . It can be tested by studying the coefficients 
E^V as we did in Sec. The rigorous bound becomes (l/2)w 2 > (1/2) (3tt 2 /8)/x max , 
which implies 

w c > 1.92 x x m 2 ax (44) 

From Sec. El we have 

u; c ~ (3.3) 2 x x~ 2 ax , (45) 

and again, the lower bound is compatible with the numerical estimate, but is not an 
accurate estimator of the actual radius of convergence. 



7. The crossover region: empirical data collapse 

In this section, we discuss the question of the interpolation between the large and 
small x m ax approximation. In Sec. 21 we have already observed that by naively using 
integral formulas derived in the large x max approximation for E^ and it was 

possible to obtain decent approximations at low x max . We will first explain how this 
works and then discuss the higher orders terms. 

In the (mathematical) limit of very small x ma x, Eq. (|23|l implies that Gq{x) ~ 
— lux 2 . Consequently, in this limit, where we do not expect the approximation to be 
accurate, the correction becomes dominant and we have 

E^(V^x max ) ~ l/iuxLj ■ (46) 
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Table 2. x (h) for k = 1 . . . 10. 



x (k) 



1 1.908 

2 2.509 

3 2.868 

4 3.152 

5 3.398 

6 3.626 

7 3.839 

8 4.042 

9 4.235 

10 4.421 



Despite the lack of justification for this formula, it is not far from the accurate answer 
derived in Sec. [SJ 

4 0) (VSJw) - 7r7(8w*L J =s 1.234/ (u;^) (47) 
Similarly for fiW, we can study Eq. 1|28[1 in the mathematical limit of small x max . 
According to Sec. tyff'(x) ~ (1 — {x / x max ) 2 ) and elementary integration yields 

E^\y^x max ) ~ (1/21) x (ux 2 max ) 2 ~ 0.0476(x^ o J 2 (48) 
Again, this is close to the accurate answer from Eq. (|3tJ|) : 

(v^Jaw) ^ 0.0411 x (w^ ox ) 2 (49) 

For higher orders, we know from Fig. El that the small Xmax approximation is not 
convergent in the crossover region. However, the shape similarities observed in Fig. 
d suggest to parametrize Rk in terms of a single function U, that can be shifted by 
a fc-dependent quantity that we denote xo(k). We have chosen xo(k) as the value of 
Xmax for which the second derivative of Rk vanishes. The numerical values are shown 
in Tableland Fig. [S] Empirically, it can be fitted quite well with 

x (k) ~0.87+1.13Vfc . (50) 

In Fig. El we show that the possibility of having 

Rk{Xmax) ^ U{x max - X (k)) , (51) 

is reasonably well satisfied for k > 2. In other words, with suitable translations, the 
Rk approximately "collapse". 

8. Conclusions 

We have shown that the calculation of the perturbative series for the anharmonic 
oscillator with a field cutoff x max can be performed reliably using numerical methods, 
and approximate analytical methods in the limit of small and large x max - For the 
coefficients of order and 1 in A, the integral formulas derived in the large x max limit , 
produce the correct leading power dependence with coefficients close to the correct 
ones, in the opposite limit (small x rna x) where they are not expected to be accurate. 
For higher orders, it is possible to approximately collapse the crossover region of the 
various order by an appropriate translation in x max . The collapse is not perfect and it 
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Figure 9. Rk(x+xo(k)) for k = 1, . . . 10. 



is necessary to resolve more accurately the details of Fig. [!5]in order to obtain accurate 
formulas for the Rk. 

We expect that it is possible to extend the small and large cutoff techniques in 
the case of higher dimensional field theory, however the only numerical method that 
can we envision for the crossover region is the Monte Carlo method (TS|. In view of 
this, it is essential to reach a proper understanding of the interpolation. 
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